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1 Introduction 

• Galaxy formation in the Universe - collapse the baryons within DMH potential wells 
(White & Rees (78)). Observational support: COBE detection (Bennett et al. (93)). 

• Formation of self - gravitating inhomogeneities of protogalactic size (Dar (95)). Origin 
of initial angular momentum (Steinmetz & Bartelmann (95)). 

• Smoothed Particle Hydrodynamics (SPH) (Monaghan (92)). TREE - SPH code (Hern- 
quist & Katz (89)). GRAPE - SPH code (Steinmetz & Muller (94, 95)). 

• Extension of our N - body/SPH method (Berczik & Kolesnik (96), Berczik & Kravchuk 
(96)). New "energetic" criteria for SF and more realistic account of returned chemical 
enriched gas fraction via SNII, SNIa and PN events. 



2 The CD-SPH code 

• The SPH code. Continuous hydrodynamic fields in SPH are described by the interpolation 
functions constructed from the known values of these functions at randomly positioned 
particles (Monaghan (92)). 
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The equations of motion for particle i are: 
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The energy equation in the particle representation has the form: 




^i(Pi,Ti) -Ai(pi,Ti) 



Pi 



The system of equations is closed by adding the equation of state: 



P i = Pi ■ (7 - 1) • «i- 



• Time integration. To solve the system of equations we use the standard algorithm of 
leapfrog integrator (Hernquist & Katz (89)). The integrator has a second order accuracy 
in the time step At. To define At we use the relation (Hiotelis & Voglis (91)): 



• The star formation algorithm. We modify the standard SPH star formation algorithm 
(Katz (92)), taking into account the presence of chaotic motions in the gaseous environ- 
ment and the time lag between initial development of suitable conditions for star formation 
and star formation itself (Berczik & Kravchuk (96)). It states that in the separate "gas" 
particle the SF can start if the absolute value of the "gas" particle gravitational energy 
exceeds the sum of its thermal energy and energy of chaotic motions: 



It seems reasonable that the chosen "gas" particle produce stars only if the above condition 
holds over the time interval exceeding its free - fall time: 



We check the number of SF acts in selected "gas" particle i. If the number of SF acts 
becomes greater than N^ x = 25 we stop any SF activity in these particles. 

We also define which "gas" particles remain cool, i.e. t coo i < tff. These conditions we 
rewrite in the manner presented in the paper by Navarro & White (93): pi > p crit . Here 
we use the value of p crit = 0.03 cm~ 3 . 

When the collapsing particle i has been defined we create the new "star" particle with 
mass m star and updated the "gas" particles rrii using these simple equations: 
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In our Galaxy on the scale of giant molecular clouds the typical values for SF efficiency 
are in the range e fx 0.01 0.4 (Duerr et al. (82), Wilking & Lada (83)). We define e as: 

1 Ef + Ef 

e = 1 — 

\Ef r \ ' 

In the code, we set the absolute maximum value of the mass of such "star" particles 
m s JZ = 10 7 M . 

At the moment of the birth, the positions and velocities of new "star" particles are 
set equal to thoes of parent "gas" particles. Subsequently, the "star" particle interact 
with the rest of "gas" and "star" particles or "dark - matter" only by gravitation. The 
gravitational smoothing length for these (Plummer - like) particles is set equal to h star . 

• The thermal SNII feed - back. For the thermal budget of the ISM, SNIIs play main role. 
Following to Katz (92) and Friedli & Benz (95), we assume that the explosion energy is 
converted totally to the thermal energy. The total energy released by SNII explosions 
(10 44 J per SNII) within "star" particles is calculated at each time step and distributed 
uniformly between the surrounding (i.e. < h star ) "gas" particles (Raiteri et al. (96)). 

• The chemical enrichment of gas. In our SF scheme, every new "star" particle represents 
a separate, gravitationally closed, star formation macro region (like a globular clusters). 
The "star" particle has its own time of birth UegSF which is set equal to the moment 
of the particle formation. After formation of these particles due to SNII, SNIa and PN 
events, return the chemically enriched gas to surrounding "gas" particles. For description 
of this process we use the aproximation proposed by Raiteri et al. (96). We concentrate 
our treatment only on the production of 16 and 56 Fe, but try to describe the full galactic 
time evolution of these elements, from the beginning up to present time (i.e. t evo \ fx 13.0 
Gyr). 

For example, if the mass of new "star" particle (with metallicity Z = 10~ 4 ) is equal to 
10 4 Mq, it produces next numbers of events: AiV SNII ~ 52.5, AAp N fx 1770, AiV SNIa fx 
8.48 during total time of evolution t evo i. 

The total masses (H, He, 56 Fe, 16 0) returned to the surrounding "gas" particles, due 
to these processes the are (in solar masses): Am H ~ 2644, Am Hc fx 881, Am Fe ~ 
8.8, Am° fx 120. 

• The cold dark matter halo. In the literature we have found some, sometimes contraversial, 
profiles of Cold Dark Matter Haloes (CDMH) in the galaxies (Burkert (95), Navarro 
(98)). For resolved structures of CDMH: phaio{f) ~ r~ lA (Moore et al. (97)). The 
structure of CDMH, as shown in high-resolution N - body simulations, can be described 
by: Phaio(r) ~ r^ 1 (Navarro et al. (96), Navarro et. al (97)). Finally, in paper by Kravtsov 
et al. (97) we find that the cores of DM dominated galaxies may have a central profiles: 
phaio(r) ~ r-°- 2 . 

In our calculations, as a first order aproximation, it is assumed that the model galaxy 
halo contains the CDMH component with Plummer - type density profiles (Douphole & 
Colin (95)). Therefore for the external force which exercises onto the "gas" and "star" 
particles by CDMH we can write as: 



Figure 1: The distribution of "star" and "gas" particles in the final step. 
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Results and discussion 

• Initial conditions. The SPH calculations were carried out for N gas = 2109 "gas" particles. 
According to Navarro & White (93) and Raiteri et al. (96), such number seems to be 
quite enough to provide qualitatively correct description of the system behaviour. Even 
such small number of "gas" particles produces a AT stor = 31631 "star" particles at the end 
of calculation. 

The value of the smoothing length hi was chosen requiring that each "gas" particle had 
Nb = 21 neighbours within 2 • h{. Minimal h m i n was set equal to 1 kpc. For "star" 
particles we use the fixed gravitational smoothing length h star = 1 kpc. 

As initial model (relevant for CDM - scenario) we took constant - density homogeneous 
triaxial configuration of gas (M gas = 10 11 M Q ) within the dark matter halo (M ha i Q = 
10 12 M ). We set A = 100 kpc, B = 75 kpc and C = 50 kpc for semiaxes of system . We 
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set the smoothing parameter of CDMH: bhaio = 25 kpc. The gas component was assumed 
to be initially cold, T = 10 4 K. 

The gas was assumed to be involved into the Hubble flow (H = 65 km/s/Mpc, ffo = 
1) and into the solid - body rotation around z - axis. We added the small random 
components of velocities (A | v |= 10 km/s) to account for the chaotic motions of 
fragments. 

The spin parameter in our simulation is A ~ 0.08. This parameters is defimed in Peebles 
(69) as: 

~ G-(M gas + M dm f^ 

If the angular momentum is acquired through the tidal torque of the surrounding matter, 
the standard spin parameter does not exceed A ~ 0.11 (Steinmetz & Bartelmann (95), 
A « 0.071°;°!)- 

• Conclusion. This simple model provides good, self - consistent picture of the process of 
galaxy formation. The dynamical and chemical evolution of modelled disk - like galaxy is 
coincident with the results of observations for our own Galaxy. Some basic distributions 
of gas and star parameters are given in figures: 

— Fig. |T[ The distribution of "star" and "gas" particles in the final step. 

— Fig. ^. V rot (r). The rotational velocity distribution of gas in the final step. 

— Fig. |3|. cr*(r), a gas {r). The column density distribution in the disks of gas and stars 
in the final step. 

— Fig. f|. T(r). The temperature distribution of gas in the final step. 

— Fig. [5] .SFR(t). The time evolution of the SFR in galaxy. 

— Fig. ||. [Fe/H](t). The age metallicity relation of the "star" particles in the "solar" 
cylinder (8 kpc < r < 10 kpc). 

— Fig. 0. A^([Fe/H]).The metallicity distribution of the "star" particles in the "solar" 
cylinder (8 kpc < r < 10 kpc). 

— Fig. g. [0/Fe]([Fe/H]). The [O/Fe] vs. [Fe/H] distribution of the "star" particles in 
the "solar" cylinder (8 kpc < r < 10 kpc). 

— Fig. § [0/H](r). The [O/H] radial distribution. 
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Figure 2: The rotational velocity distribution of gas in the final step. 
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Figure 3: The column density distribution in the disks of gas and stars in the final step. 
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Figure 4: The temperature distribution of gas in the final step. 
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Figure 5: The time evolution of the SFR in galaxy. 
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Figure 7: The metallicity distribution of the "star" particles in the "solar" cylinder (8 kpc 
< r < 10 kpc). 
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Figure 8: The [O/Fe] vs. [Fe/H] distribution of the "star" particles in the "solar" cylinder (8 
kpc < r < 10 kpc). 
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Figure 9: The [O/H] radial distribution. 



